Spin glass transitions of smectic-A crosslinked elastomers 
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Elastomers are artificial polymeric materials created for industrial and commercial applications. 
Depending on their purpose, they are performing in different species and structure modifications. 
Our studies focus on the systems of elastomers randomly standing-distributed in a smectic A (SmA) 
liquid crystal. Basing on the suggestion following from the experiment, that at a phase transition 
from Smj4 to nematic phase caused by an increase of a crosslink concentration, such a system 
survives a percolation transition at low crossilink concentrations, we propose a modeling explaining 
this phenomena. We approve the three-dimensional Villain spin glass model and apply lattice Monte 
Carlo (MC) techniques on differential forms on a dual lattice, that is an alternative of a replica trick, 
developed for nematic elastomers in the 3D XY universality. In the results of that we have confirmed 
a concentration phase transition of percolation nature at a small crosslink concentration (~ 10 weight 
%). 
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INTRODUCTION 

A great number of quantum phase transitions in poly- 
mers and liquid crystals exhibits a percolation phe- 
nomenon. 

Particularly, percolation properties of the transi- 
tion from smectic-^ (SmA) to nematic [N) phase in 
crosslinked smectic elastomers were observed on the ex- 
periments on X-ray [1] and ^H-NMR spectroscopy |3|. 
Percolation in such systems may occur at the increase of 
the crosslink density and/or at the increase of tempera- 
ture as well. 

Elastomers are polymer networks composed by 
crosslinked polymer chains. They consist of nodes and 
links obeying ascertained rules of self-organization [^ 
|6|. Polymeric random network systems are differing 
with types of ordering, orientational and translational, 
caused by topological defects (dislocations, vortexes etc.) 

For mere polymers, the percolation problem is often 
discussed in the light of the analogy to dimers and con- 



nective conductors, proposed by de Gennes [9|, [10|, [llj, 
and fractal dimensionality. However, crosslinks are not 
interacting via only central forces (refs. in [1^) (which 
correspond to a bending energy), but both strength and 
bending energies are presented, carrying out to a high 
percolation threshold. At any modifications in a poly- 
meric network, using of the de Gennes's theory may be- 
come more embarrassing. 

In mere smectics, the linear nontrivial disclinations 
(3D vortexes or dislocation loops) as well as point 
monopole-type defects take place 13l-ll5|. From a rep- 
resentation of topological defects in lattice theories, per- 
colation associates with condensation of vortexes in the 
high temperature deconfinement phase [l6|. Where, sim- 
ilarly that in polymers, the fractal-dimension domains of 



this phase will exist on a non-integer lattice, the Polyakov 
line plays the role of the order parameter. 

For pure smectic-^ liquid crystals, moving by the sec- 
ond order tilt transition to the nematic phase, Dasgupta 
proposed the loop inverted analog of the superconductive 
XY model [15j . The action of this model is equivalent of 
the Villain's action |17|, if in the anisotropic XY model, 
the de Gennes coefficient Ki equals to zero 15|, llSl, ll9| . 
Dasgupta ([l5| and references therein) derived mapping 
of the partition function of the lattice de Gennes model 
onto that of a dislocation loop model. A vortex, as well as 
a dislocation with its similar properties, expresses a non- 
trivial classical minimum of the action [13| of a scalar 
field theory $ with broken global U{\) symmetry, and 
the XY model is its quantized variant. 

Relating to the 3D XY universality, percolation tran- 
sitions could be described via the U{\) gauge field A [20| 
on a dual cubic lattice at critical temperature /3c ~ 0.4542 
in appointed fractal dimensionality I?/, while quantita- 
tive characterization of the order parameter is remaining 
complicated by ambiguity in the field A [l5|, [2ll |. 

In a general case, for spin glass models of polymer net- 
works, one applies the replica methods [22l - l25| . 

With understanding the SmA-A^ transition in alkyl- 

cyanobiphenyl nCB (in particular, in 8CB) caused by 

porosity, it has become possible to subdivide the theory 

of Sm^-A^ transitions into the three schemes (Anopore, 

I — I I — li 
silica Aerogels and Vycor glass ) [14 I26l | . As was shown 

numerically, in bulk nematics [26[, due to the coupling 
between smectic and nematic order parameters and direc- 
tor's fluctuations, the SmA-A^ transition deviates from 
3D XY universality to a nonuniversal crossover value be- 
tween 3D XY and a tricritical behavior with the critical 
exponent a of the effective heat capacity {a ranging be- 
tween 0.26 and 0.31 confinements). 

In the lack of elastomers, in the 3D superconductive 



universality [ij, |27| with a scalar gauge field A: = 0, a 
dimensionality of defects (vortices) isD — k — 2 = 1 
|20| . and may be successfully assigned to aerogels with 
the smectic ordering. A further fractal contribution into 
the additional material (defectless) dimensionality may 
be produced by the crosslinked elastomers, confined in 
such smectics. 

So the order parameters of both the smectic-^ and the 
polymer network is to be conflated as a whole. In this 
paper, we present the SmA-N transition in the system of 
smecticA and crosslinked elastomer as behaving in frames 
of the Villain model and specify the spin-glass treatment 
of it. 



THE VILLAIN MODEL ON A DUAL LATTICE 

Except the continual theory of the thermotropic SmA- 
N transition in crosslinked elastomers [29|, the related 
thermod yna mical expression in similar SmC systems was 
given in [30j. 

Referring to the continual model [3l|, the authors of 
[3 01 have found, that though in lack of the external mag- 
netic field, the expressions of the free energy density for 
the smectic and nematic elastomer are analogous, how- 
ever, the nematic elastomers possess so called "soft elas- 
ticity" , e. g. there is no homogeneity of a nematic direc- 
tor field, but splay, bend, and splay-twist fluctuations are 
presenting together. This fact, as well as the conclusions 
of [1^ from the " spring-ball" numerical modeling of gels 
and filled polymers, being crowned with a non-affinity in 
the deformations here, compels us to accept a proposi- 
tion of [31[, that spin glass models are only satisfying to 
receive statistical values revealing the content dependent 
percolation threshold in the SmA-crosslinked elastomer. 

Then we will introduce an assignment of lattice spin 
variables to any components and configurations of this 
system. On its Abelian discrete Zff group exp{2jfc7r/n}, 
k = 0,...,n — 1, an oriented link bordering plaquettes 
with the appointed orientation and non-orientable sur- 
faces are permitted for even n. 

Vortices of the Villain model are decoupled with spin 
waves, but the XY model has the same topological char- 
acteristics, as the Villain model. It is reasonable to as- 
sociate frustration effects with a crosslink configuration, 
i.e. with unfrustration networks, where the thermal fiuc- 
tuation does not destroy the long range order, according 
to the spin glass topology in three and more dimensions 
32l l33j . As to the percolation phenomena in these elas- 
tomers at the SmA-N transition, a random field replaces 
onto a boundary effect on the unfrustration network. 

Further, we follow these ideas on a frustration for spin 
glasses [32|, |33| : " a frustration network is dual to a set of 
the (d — 2)-cells in the dual lattice which form {d — 2)- 
dimensional complexes which are closed or whose bound- 
aries end at the boundaries of the dual lattice A* to the 



lattice A." 

From the expression 



n03p= n (4.)'=i, 



(1) 
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one follows that an even number of the plaquettes must 
be frustrated in a cube c, and the complexes of dual 
{d — 2)-cells to the frustrated plaquettes can not end in 
any cube. Here in ([T]), J^j = Ji,j/\Ji,j\ mean the signs 
of the couplings Jij in the Toulouse's introduction of 
frustration for a loop c along the bonds of the lattice 

('/'(c) ^U{t,j)CcJ^,j)- 

Percolation in such cases is formulated [32] as a 
Bernoulli bond percolation process on the Z^ lattice with 
densities x and (1 — x), respectively for the sorts of spins, 
the bonds of negative couplings Jij percolate for x near 
1/2 (just like in Z'^). The thermal fluctuation does not 
destroy the long range order in the unfrustration net- 
work. Namely we can expect the existence of the spin 
glass phase on the Z^ lattice for the density x near 1/2. 

In |32| . there are considered the cube configurations 
whose two-cells are made of either some unfrustrated 
pairs of frustrated and some frustrated plaquettes, or all 
made of unfrustrated pairs, so that the way of choosing 
unfrustrated pairs as two-cells is not unique, and in Z^ 
the ground state is highly degenerate; one may only hold 
a condition that any loop which is made of the boundary 
bonds of the unfrustrated pairs should be homologous to 
zero in the cube. 

Related studies for all frustrated plaquettes has been 
numerically provided in Z^ for a large dice lattice 
(L =84) \m\. 

The Hamiltonian of the XY spin glass model on a cubic 
lattice is written as follows: 
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This model supports vortices 35| unbind at high temper- 
ature via the Berezinskii-Kosterlitz-Thoules (BKT) tran- 
sition. Here 9i indicates the angle of the XY spin at the 
i-the site and the summation is taken over all nearest 
neighbor pairs, and the coupling constants Jij obey such 
a distribution as the ± J or Gaussian function as in the 
Ising spin glass. We set the Boltzmann constant fc^ to 
unity. 

On the lattice bonds, Aij = —Aji, the magnetic cur- 
rents connect the centers of lattice boxes and obey 



^'' " To 



drA(r). 



(3) 



These currents do not fluctuate, their sum belong the full 
lattice perimeter equals to 27r/, J2kvadratic ^v = ^ti"/, 



the equality / 



00 



expresses the magnetic quanta 



and the magnetic current per one lattice cell, in fact, 

/e[0,i]. 
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FIG. 1: Temperature dependence of an average energy at 
different crosslink concentrations 



FIG. 2: Temperature dependence of an average heat capacity 
at different crosslink concentrations 



The substitution e'^™"'^' 



-> 



E- 



-C(ipi3-27rri)^/2 



f[35j. where ipij = 9j—9i is the rotation angle of neighbor- 
ing spins on a lattice, and n is whole) made by Berezin- 
skii and Villain (see references in [36], and [37]) allows 
us to write the spin glass partition function of ^ on the 
original lattice: 



ve 



E 



exp{[- 



J, 



-ei + 27rp,,f]}, (4) 



here pij are winding numbers lengthwise a loop, 6i and 
^j mean phases in given lattice sites. Vd denotes the 
integral over all link variables 6''s. 

Possessing view of the Coulomb gas partition function, 
(J4|) corresponds with the partition function derived in [l9| 
on a dual lattice describing a set of integrating vortex 
loops defined in the generalized Villain's model with a 
non-zero " chemical potential" . 

With the technique of differential forms 



la 33 on the 



lattice, one builds the dual lattice configurations corre- 
sponding to ([4]) with the next partition function 

zduai^ ^ j2 iexp{-A'||dfc + ;||2-/?'||;||2}, 

(5) 
here the field of "currents" I attached to the links and 
flowing along the C^ circuit, and the "monopole" field 
k attached to the sites, which is provided with currents 
along C° . The expression ([5]) also has an equivalent sense 



as ^ r 



Z„ 



'■^compact ^defects ' ^ jionconipact • 

Further connection between original and dual lattice 
variables are performed according to the transformation 
[l6| in the gauge invariant form of Z^ space. The terms 
Jexp{-A'||dA: + Zf - /3'\\l\\^} reflect the exp{-4^(6'j - 
9i + 2TTpij)^} 's sense after the differential form trans- 
formation. Coulomb interactions between currents and 
monopoles are straight to the linking number of currents 
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FIG. 3: Configuration dependence of probability from the MC 
algorithm on the dual lattice at different crosslink concentra- 
tions, 0.06, 0.10, 0.15, 0.20 % 



NUMERICAL RESULTS AND DISCUSSION 

In our MC modeling, we choose a lattice of size 48'^ and 
three type of whole spin variables, corresponding to smec- 
tic molecules (0), polymer networks (1) and crosslinks 
(-1), all occupied a cubic lattice sites. Closed Wilson 
loops and open loops appear in a system. Well-known 
standard loop algorithms were involved ^], ^j. The 
specific heat expression was improved in accordance to 
the Villain action properties [41] (Figs 1-3 confirm, that 
temperature dependent transitions are not percolating). 

For the four-dimensional Z^ gauge theory, there was 
obtained the transformation connected the parameters 
/3', A' and /z' of the dual lattice with the original ones, 
which is follows [16|. 

In the 3D case (/i', /i equal to zero), the reverse trans- 
formation for these coefficients of (5) is performing by 
the equality (3 = -^ and vice versa. 

Fig. 4 illustrates extremums of an average heat ca- 
pacity, which appear when a crosslink concentration sur- 
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FIG. 4: Average heat capacity versus a crosslink concentra- 
tion on the dual lattice at different temperatures in units 
/3 = ^ (inset). Error bars are not shown. For visualization of 
a heat capacity extremums, a Lorentzian fit for /3 = 0.2021 is 
shown 
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FIG. 6: The correlation function versus temperature on the 
dual lattice at different temperatures. Error bars are not 
shown. 
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FIG. 5: Average susceptibility versus temperature on the dual 
lattice at different temperatures. Error bars are not shown. 



FIG. 7: The correlation function versus crosslink concentra- 
tion on the dual lattice at different temperatures. Error bars 
are not shown. 



passes 10 %. This result confirnis the experimental data 
of [Ij , but does not satisfy a confirmation on percolation 
of such a transition: the configuration averages and their 
probabilities shows (Fig. 3), that they meet at tempera- 
tures close to zero, but not at /3 =0.4542 or 0.5, predicted 
in [20j an 32 1 in frames of the gauge BKT and spin glass 
theories. 

It is known, that long-range ordering may be described 
in terms of spin glass models 3]| , 42| by means of the 
replica method, but, whereas our model is of nearest- 
neighbor interactions, we used a loop model. In 31j . 



there was predicted a phase transition in nematic elas- 
tomers via an order parameters, expressed as a magne- 
tization, but in the version of a model as ([5]), this value 
does not exist. 

Roughly speaking, using a monopole part of the 
Abelian Polyakov loop, we have got also correlators 
(Fig. 7) due to the standard procedure [43|, |4J] as func- 
tions of a crosslink concentration, which do not directly 



confirm percolation, but manifest features increasing at 
growing temperatures at the 10-percent concentration of 
crosslinks. 



CONCLUSION 

To resolve the problem of percolation properties of 
elastomer configurations, we applied the technique lat- 
tice differential forms, builded on the Hamiltonian of the 
3D XY spin glass model, the main advantage of which is 
reduction of a dimensionality of the system and making 
an ability to estimate its thermodynamic values without 
cumbersome calculations. This method is alternative to 
the replica approach. However we have confirmed the ex- 
perimentally observable percolation concentration phase 
transition at an increase of a crosslink concentration ex- 
ceeding 10 percents. 
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